Curvature in conformal mappings of 2D lattices and foam structure 
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The elegant properties of conformal mappings, when applied to two dimensional (2D) lattices, 
find interesting applications in 2D foams and other cellular or close packed structures. In particular 
the 2D honeycomb (whose dual is the triangular lattice) may be transformed into various conformal 
patterns, which compare approximately to experimentally realisable 2D foams. We review and 
extend the mathematical analysis of such transformations, with several illustrative examples, and 
an account is given of the related work in energy minimisation problems. New results are adduced 
for the local curvature generated by the transformation. 



I. INTRODUCTION 

The relationship w = f(z), where / is any analyt- 
ical function, can be viewed as a mapping which sets 
up a correspondence between the points of the z and w 
planes. Such mappings are known as conformal map- 
pings. The geometrical operations of inversion, reflec- 
tion, translation and magnification are all examples of 
conformal transformations in Euclidean space. Confor- 
mal mappings have a number of interesting properties, 
the most important being isogonality: any two curves 
that intersect are transformed into curves that intersect 
at the same angle. 

Consider a discrete 2D set of points in the z-plane gen- 
erated by two primitive vectors: the resulting structure 
in the image plane, due to a conformal mapping, is known 
as a conformal lattice. It is therefore a purely geomet- 
rical object. A strictly conformal crystal is a physical 
system consisting of particles located on the sites of a 
conformal lattice. A conformal crystal is a physical sys- 
tem, in which the arrangement of particles approximates 
a conformal lattice 

There are numerous examples of conformal crystals oc- 
curring both in nature and in the laboratory; despite 
this the geometric properties of conformal crystals are at 
present poorly understood. In the present paper we ask: 
what determines the local curvature in these conformal 
patterns? We shall begin our analysis by considering the 
equation for the complex curvature, which was derived by 
by Needham Q and also by Mancini and Oguey (1J fjj. 
For a given line in the z-plane, we shall see that whereas 
the first derivative of the transforming function relates 
the direction of a line to its transformed counterpart, the 
induced curvature involves the second derivative. Here 
we shall use the equation for the induced curvature to 
compute the mean and mean square curvature of the con- 
formal lattice 

One of the most easily recognised examples of a nat- 
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urally occurring conformal crystal is the phyllotactic de- 
sign of a sunflower Q . Another example is the so called 
"gravity's rainbow" structure, which is the name given 
to the striking arrangement of arches formed by a clus- 
ter of magnetised steel balls in an external force field 
More recently, conformal lattices have been shown 
to have a connection with disclinations in 2D crystalline 
structures ® 0- 

Conformal crystals can be physically realised by sand- 
wiching an ordered, quasi 2D, foam in a Hele-Shaw cell 
with non-parallel plates, see for example [8J. Yet another 
method involves the use of ferrofluid foams in magnetic 
fields @. The advantage of these foam-based methods 
is that a variety of conformal crystals can be realised by 
tuning the geometry of the experiment. However, the 
use of foams to approximate conformal lattices involves 
two complications which have to be taken into account. 
Firstly, as was shown by Mancini and Oguey (see [3| [J]) 
the curvature of the soap films perpendicular to the glass 
plates has to be taken into account (in section III we 
summarise the main results of the work conducted by 
Mancini and Oguey). Secondly, we compute the higher 
order corrections in the expression for the complex cur- 
vature and show that is only when these terms are small 
that the foam is a good approximation to a conformal 
lattice. 

The paper is organised as follows. In section II we in- 
troduce the complex curvature and give some properties 
of conformal transformations. The relationship between 
ordered 2D soap froths and conformal transformations 
is detailed in section III. In section IV we calculate the 
mean curvature and the mean square curvature of the 
conformal lattice when the original lattice in the z-planc 
is free of curvature, for a purely 2D system. We illus- 
trate these results with some examples which include the 
case of complex inversion. In section V we generalise our 
results and include the case where the original lattice in 
the z-plane has a curvature. In section VI we compute 
the higher order terms in the expression for the complex 
curvature. In section VII we discuss some further appli- 
cations for conformal crystals, these include disclinations 
in 2D crystals and the intelligent design of 2D computa- 
tional meshes. Finally, we give our conclusion in section 
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VIII. 



II. 



SOME PROPERTIES OF CONFORMAL 
TRANSFORMATIONS 

A. Scaling of Areas 



Although a conformal mapping w = f(z) preserves 
angles (the isogonal property) it does not preserve areas. 
If ds z — (dx 2 + dy 2 )? is a small element of line in the 
(x,y) plane, upon being mapped to the w-plane it will be 
magnified and have a length given by, 



ds v 



dw 
dz 



ds z = |ds 



(1) 



where f n (z) is the nth derivative of the function f(z). 
Hence a small element of area in the z plane, denoted 
by dA z , will upon being mapped to the w-plane have an 
area 



dA„ 



dw 
dz 



dA z = \fi(z)\ 2 dA z 



(2) 



B. Complex Curvature 

Consider a curve K in the z-plane: if we apply an ana- 
lytical mapping f to this curve then it will transform into 
another curve in the image plane, which we denote by K. 

Let us now choose some arbitrary point on K which 
we denote by p = x + iy = re 1 ® . The unit tangent vector 
to the curve at point p is given by £(<f>) = e 1 ^, where <fi 
is the angle the tangent vector makes with the x axis; 
for an illustration see the top part of Fig. (JTJ). Upon 
applying an analytical transformation, point p is mapped 
to a new point in the image plane which has coordinates 
f(p) — u + iv — Re iA . It has been shown, see 0, that if 
the instantaneous curvature of K at the point p is given 
by re, then the instantaneous curvature at f(p) is given 
by, 



1 



\h(p)\ 



Im 



MP) 



\A(p)\ 



(3) 



where Im is the imaginary component, see the bottom 
part of Fig. (JTJ) . The first term in Eq. {3]) is the curvature 
in the image plane if K is a straight line. If however K 
has a curvature re, then in the image plane this additional 
curvature is scaled by a factor of l/\fi(p)\; again see @ 
for a beautiful derivation of these and other results. 



Example: f{z) = z a 

The utility of Eq. ((3]) can be demonstrated by a short 
example. Consider the mapping, 



w 



(4) 
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FIG. 1: Top: For simplicity we choose to consider a straight 
line K, in the z-plane, which has no curvature (i.e. k = 0). As 

shown the tangent vector £((/>) at the point p makes an angle 
of (j) with the x axis. Bottom: In the image plane the effect 
of the analytical mapping is toyield a curve K, the curvature 
of which is given by k. After 



where 



/i(z) = az a ^ and f 2 (z) = a(a - l)z Q ~ 2 , (5) 



and let us represent the z-plane and the image plane in 
terms of complex polar coordinates so that, 



and w — Re 



Upon substituting Eq. (O into Eq. ([3]), we find that the 
the effect of mapping a straight line in the z-plane (i.e. 
re = 0) is to yield a curve in the w-plane with curvature, 
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ar u 



a — 1 
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Writing out explicitly and simplifying, this becomes, 
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yOL — 1 



Im 
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Im 
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sin(<£-0). 



Note, the curvature is expressed in terms of the coordi- 
nates of the z-plane. To get the curvature in the w-plane 
we must use the relationship R — r a , which gives 

~ a - 1 1 . , , 



a R v 

Thus the maximum curvature occurs when the tangent 
vector is perpendicular to the vector connecting the ori- 
gin to the point p and is given by, 



1 1 



R 



(6) 



This means that lines, in the z-plane, drawn perpendic- 
ular to the vector connecting the origin with point p will 
acquire the greatest curvature, while lines drawn parallel 
to it will not suffer any curvature. 



III. PROPERTIES OF 2D SOAP FOAMS AND 
THEIR RELATIONSHIP TO CONFORMAL 
TRANSFORMATIONS 

A dry 2D foam consists of 2D bubbles separated by 
lines which meet at vertices. Such a foam can be repre- 
sented by a network consisting of 2D cells separated by 
ID edges. Equilibrium conditions impose strict restric- 
tions on the topology and geometry of such a foam net- 
work (see 10]): Plateau's law stipulates that the edges 
can intersect only three at a time and must do so at 
an angle of 27r/3; the edges themselves have a constant 
curvature, that is they form circular arcs, and the curva- 
ture of the arcs is related to the corresponding pressure 
difference between the adjacent bubbles by the Laplace- 
Young relation. As a consequence of the Laplace- Young 
relation it follows that the sum of the curvatures of the 
three edges at a given vertex vanishes, or equivalently 
that the mean curvature of the adjoining edges vanishes. 

All of the above conditions are automatically satisfied 
by complex inversion, f(z) = 1/z, or more generally a 
bilinear conformal transformation. A bilinear transfor- 
mation has the form, 



az 



It can be decomposed into four sequential transforma- 
tions: translation, inversion, expansion and rotation, and 
a final translation, for details see [2| ■ Of these four steps 
only inversion is not trivial. Conformality ensures that 



the mean curvature at each vertex vanishes, while only 
complex inversion (or a bilinear transformation) has the 
special property that it will map a circular arc into an- 
other circular arc. Given a dry 2D foam structure at 
equilibrium, inversion will therefore produce a new equi- 
librium structure; this was discussed by Weaire who 
used it to provide a neat proof of the Decoration Theo- 
rem. 

What we refer to as a quasi 2D foam is a foam consist- 
ing of a single layer of bubbles trapped between a pair 
of narrowly separated, parallel, glass plates [ll| - such a 
confining geometry is known as a Hele-Shaw cell. Pro- 
vided the separation of the glass plates is below a critical 
value the quasi 2D foam is a good approximation to a 
2D foam, i.e. the equilibrium conditions for 2D foams 
as stated above - Plateau's law and the Laplace- Young 
relation - also apply to quasi 2D foams to a good approx- 
imation. If the separation of the plates is greater than 
this critical thickness the bubbles undergo a three dimen- 
sional instability and the foam is no longer a monolayer 

It is important to remember that the boundaries be- 
tween bubbles in such quasi 2D foams are indeed 2D soap 
films and not ID edges as they are often approximated. 

In the ideal case, if all the bubbles trapped in the Hele- 
Shaw cell have the same volume (monodisperse) , then an 
ordered quasi 2D foam can be realised. Upon viewing 
the Hele-Shaw cell from above (i.e. from the direction 
perpendicular to the plates) the bubbles are observed to 
form a honeycomb structure. 

Let us now assume that the bottom plate of the Hele- 
Shaw cell is flat while the upper plate is slightly an- 
gled or curved, if the bubbles are monodisperse, then 
this imposes a specific variation in the area of the bub- 
bles. The variation in the bubble area can be made to 
closely match the variation required by a given confor- 
mal transformation - as stipulated by Eq. (5]). Note that 
although the area of the bubbles (as observed from the 
direction perpendicular to the bottom plate) may change 
their volume remains constant, thus the height of the up- 
per surface is related to the analytical function f(z) by 
h(w) — l/dA w ). This fact has been used to transform 
the (straight-edged) honeycomb structure, in a variety of 
ways, to generate approximations of conformal lattices 
Q (i.e. conformal crystals), see for example Fig. @. 

However, the use of quasi 2D foams to simulate con- 
formal lattices leads to a contradiction: on the one hand 
equilibrium conditions demand that the curvature of the 
soap films is constant; on the other hand only bilinear 
transformations have the property that the transformed 
edges have a constant curvature. Thus, with the excep- 
tion of the bilinear transformation, the resemblance of 
any quasi 2D foam to its corresponding complex map 
may be questioned. The resolution to this problem was 
recently provided by Mancini and Oguey, see Q [1], by 
considering the curvature of the soap films perpendicular 
to the plates of the Hele-Shaw cell. 

Here we briefly state the main results derived by 



FIG. 2: From a realisation of the logarithmic map (also known as the gravity's rainbow structure), (a) Experimental setup. 
Monodisperse bubbles are trapped between two glass plates, the lower plate is horizontal while the upper plate is angled. This 
arrangement induces specific variation of the bubble area with position in the foam, (b) Numerical mapping of the perfect 
honeycomb using the transformation w — (ia)' 1 log(iaz). (c) Experimentally obtained pattern. The shading shows how lines 
that were initially straight are transformed. 



Mancini and Oguey; we refer the reader to their pub- 
lications for a full account. Fig ([3]) shows a soap film 
(dotted line) in a Hele-Shaw cell with non-parallel plates. 
The height of the top plate, with respect to the bottom 
plate, is described by a function h(x\, xz)- 

The top part of Fig ([3]) shows the cell from above (i.e. 
from the direction that is perpendicular to the bottom 
plate). It can be seen that the dotted line, which rep- 
resents the soap film, has a curvature in the direction 
parallel to the bottom plate. We shall call this the longi- 
tudinal curvature and its instantaneous value at the point 
q on the curve is denoted by At. Fig also shows the 
normal n and tangent £ unit vectors to the curve at the 
point q. 

Now consider the section of the soap film at the point 
q as observed from the direction parallel to the bottom 
plate - this is shown in the bottom part Fig ^ - where 
the upper surface can be assumed to be locally straight by 
a sufficient magnification close to the point q. As recog- 
nised by Mancini and Oguey if a soap film is sandwiched 
between two non-parallel plates then, by the normal in- 
cidence condition, the film must bend in the direction 
normal to the plates. We shall call this curvature, which 
is in the direction perpendicular to the bottom plate, the 
transverse curvature. By assuming that both h and |V/i| 
are small, where Vh locally indicates the direction of 
maximum slope, to lowest order the transverse curvature 



of the soap film is can be shown to be given by 

_ Vft 

K t ~ -11- — 

h 

Thus Laplace's law can be stated as 

P\ — P% = 2jH = 7(k + Kt) (7) 

where P\ and P2 are the internal pressures of the two 
bubbles either side of the soap film. Here H is the total 
curvature of the soap film. 

Mancini and Oguey considered two cases: the first is 
the trivial case when the pressure is the same for all the 
bubbles in the conformal crystal (although this case is 
somewhat artificial from an experimental perspective), 
the second is where the volume is the same for all bubbles. 

For the case of constant pressure it is readily apparent 
that Kt = —At, thus the transverse and longitudinal cur- 
vatures cancel and H = 0, this is true for every soap 
film and at every point on the soap film. 

If on the other hand the volume is constant then it is 
found that (see and [J) 

Kt = —Ik. 

The the total curvature, of the soap film, is now not 
constant but is given by 

H = -(« + K t ) = ~2 K - 
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FIG. 3: A soap film (dotted line) confined between a pair 
of non-parallel glass plates. Top This shows the soap film 
as observed from the direction perpendicular to the bottom 
plate. For an arbitrary point on the curve, labelled q, we 
can define the normal n and tangent £ unit vectors. The 
transverse curvature at this point is denoted by k. Bottom 
After [| and @|. This shows the soap film as observed from 
a direction parallel to the bottom plate. Since the film meets 
both the top and bottom faces orthogonally it must curve, to 
lowest order this curvature is given by n t ~ —n ■ Vh/h 



This result is remarkable in that for the constant volume 
case, provided that both h and |V/i| are small, the to- 
tal curvature is related in a simple way to the transverse 
curvature - which can be readily observed in an experi- 
mental situation. Again this result holds for every soap 
film and at every point on the soap film. Laplace's law, 
given by Eq. ([7]), now reads 



IV. LATTICE CURVATURE AND CONFORMAL 
TRANSFORMATIONS 



For the remainder of this paper we shall only be con- 
cerned with complex curvature as defined by Eq. ^ . In 
terms of a quasi 2D soap foam confined in a Hele-Shaw 
cell with non-parallel plates this is the longitudinal cur- 
vature, as defined in the previous section. 

In this section we shall examine the curvature that a 
lattice, in the z-plane, acquires upon being mapped to 
the w-plane. Although the triangular lattice is of pri- 
mary importance, the results derived are general enough 
to include other structures. This includes the square lat- 
tice and the honeycomb structure. 

Consider a vertex in the z-plane located at p=x+iy: 
connected to this vertex there are, n = 0, 1,2. ..to > 3, 
straight edges. The angular separation between the tan- 
gent vectors of successive edges being equal to 27t/to. See 
Fig. (fj|. In the case of n = 6 or n = 4 it is possible to 
define two edges, which are connected to the vertex, as 
primitive vectors and thus tessellate the z-plane forming 
triangular or square Bravais lattice, respectively. 

We cannot, however, generate the honeycomb struc- 
ture in the same way (using a vertex with n = 3) since 
the honeycomb structure is not a Bravais lattice. Nev- 
ertheless we can tessellate the z-plane with a honeycomb 
structure, in a unique way, by making a Voronoi con- 
struction of the points which define a triangular lattice. 
This is because the triangular lattice and the honeycomb 
structure together define a Voronoi/Dclaunay dual, as 
shown in Fig. ([5]). 

In each case the tessellation in the z-plane can be 
thought of as consisting a vertex to which their are con- 
nected a number of straight edges of finite length (6 edges 
for a triangular lattice, 4 edges for a square lattice and 
3 edges for a honeycomb structure). In the following 
we shall examine the effect a conformal mapping has on 
these edges. We shall compute their curvature in the im- 
age plane - from which we can compute the mean and 
mean square curvature of the edges connected to a given 
vertex in the image plane. 

To keep things simple we only consider the case where 
the original edges in the z-plane, connected to a given 
vertex, are free of curvature. This condition will be re- 
laxed in the next section. 



Pi - P 2 = 2 1 H 



Since the total curvature is not constant it is important 
to know by how much the mean curvature deviates from 
being a constant (as required by the conditions of equilib- 
rium). This discrepancy sets a limit on the applicability 
of conformal transformations in the context of foams. In 
section VI we calculate the magnitude of this discrepancy 
to lowest order. 



A. The Complex Curvature of an Edge Upon 
Being Mapped to the Image Plane 

Upon applying an analytical mapping, the vertex will 
now be located at a point f(p)=u+iv in the w-plane. 
Since all the edges are straight (i.e. K n = 0), using Eq. 
([3]) it is found that the m edges have a new curvature 
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given by, 



Im 




Re 



FIG. 4: Consider a vertex located at point p in the z-plane, 
connected to which there are m straight edges. From the 
point p we can draw tangent vectors along each edge which 
is denoted by £o,£i-..£m-i. The angular separation between 
the tangent vectors is given by 2ir/m. We define an angle 7 
which is the angle between the x-axis and the first vector fo- 
Thus the angle that the nth vector makes with the x axis is 
2im/m + 7. Clearly the edges are invariant under a rotation 
of 2n/m and so we restrict 7 to the range < 7 < 2-k jm. 




FIG. 5: The white circular points in this diagram represent 
the lattice sites of a (partial) triangular lattice. The (par- 
tial) Voronoi construction partitions the plane into cells each 
enclosing a single lattice site, where each point inside this 
cell is closer to that lattice site than to any other. In this 
diagram the cells generated by the Voronoi construction are 
drawn using heavy lines, it can be seen that the Voronoi cells 
form a honeycomb structure. The honeycomb structure can 
be thought of consisting of a number of vertices (black circu- 
lar points) to which are connected three edges. Also shown in 
this diagram is the Delaunay triangulation (drawn using light 
lines) for the same lattice sites. The Delaunay triangulation 
is defined as the geometric dual of the Voronoi construction. 
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where 



and 



\fi(p)\ 2 

/ 2 (P) 



(8) 



T 



hip) 



= 02-01 = Arg 



/ 2 (P) 



MP) 



and Arg stands for the complex argument. From here on 
we shall adopt the notation that 

6n = Arg[f n (p)]. 

The nth tangent vector is given by, 

f - e ^+i) 

sn c i 

where 7 is the angle made by the Oth tangent vector, in 
the z-plane, with the x axis. This is shown schematically 
in Fig. ([4]). Eq. (|SJ) can be written as, 



\hip)\ . f n 2™ 

■ sin I fc> H h 7 



\hipW 



\hip)\ . , 

■ sin s 



\MpW 



m 
m J 



(9) 



where s = 7 + O . 



B. The Mean Complex Curvature of the Edges, 
Connected to a Given Vertex, in the Image Plane 

Let us define the mean curvature of the edges in the 
z-plane as, 



^ n—m—l 
m — 1 



n=0 



similarly, the mean curvature of the edges in the w-plane 



^ n=7n— 1 

Cm / ■ 

m ' ' 



Substituting Eq. © into Eq. |(TUJ) yields, 

n=m — l 



1 1/2Q01 

TO |/i ip)\ 2 



£ sii 



n=0 



2irn 



(10) 



(11) 
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Using the theorem, 



First Example: 



,1/2 



n—rn— 1 

E 

n=0 



sin 



2irn 
m 



= 0, 



(12) 



it is found that the mean curvature of the edges, con- 
nected to a given vertex in the image plane, vanishes. 
This result is rather trivial it merely expresses the fact 
that in the image plane, connected to each vertex, there 
are pairs of opposed edges for which the sum of their 
curvatures always vanish. 



C. The Mean Square Value of the Complex 
Curvature of the Edges, Connected to a Given 
Vertex, in the Image Plane 

Since the mean curvature always vanishes, the mean 
square curvature of the edges in the z-plane is simply 
defined as, 



-. n—m — l 

- y k 



2 

n i 



n=0 



and in the w-plane the mean square curvature is, 

n— m— 1 



— 2, kr, 

m * — ' 



n=0 



Substituting Eq. © into Eq. dH} yields, 

2 n—m — l 



rr 1 ( 1-^)1 



m \ \h(p)\ 2 
Using the theorem, 

n—m — l 

E * 

n=0 

Eq. P^|) simplifies to give, 



E. o / 27m 
sm 2 s H 
V m 

n=0 v 



sin 2 f s + | 



27T7T, 
TO J 



m 



^_1/ I/2(P)I 

Vm ~2\\f 1 (p)\* 



(13) 



(14) 



(15) 



(16) 



(17) 



Thus the mean square curvature of the edges, connected 
to a given vertex in the image plane, is independent with 
respect to the orientation of the vertex in the image or z 
planes. 



Examples 

To illustrate the results derived for the mean and mean 
square curvature, we now examine two conformal lattices 
which have been generated by applying a transformation 
to a circular disk cut out of a triangular lattice. Just such 
a cut out is shown in the top left hand part of Fig. ([6]), 
it can be seen that the lattice is free of curvature. We 
examine the case of w = z 1 / 2 and complex inversion w = 
l/z. 



Applying the mapping, 



= ,1/2 



(18) 



to the triangular lattice in the z-plane yields a conformal 
lattice in the image plane, see the top right hand side 
of Fig.©. Note, the mapping has been produced by 
allowing the range of the z-plane to extend to A < 4ir. 

To find the mean square curvature we substitute Eq. 
(fI8"]l into Eq. ([I7|). which yields 



(19) 



where we have transformed from the coordinates in the 
z-plane (r, A) to the w-plane coordinates (R, A). 

Second Example: l/z 

Similarly, the effect of inversion is shown in the bottom 
part of Fig.©. To compute the mean square curvature 
substitute 

f(z) = l/z. 
into Eq. (Tl"T|) , which yields, 



r 
~2 



2 



(20) 



Comparison with numerical results 

To verify our results we estimate the mean square cur- 
vature directly from the transformed lattice. To do so, 
note that any given hexagonal cell in the image plane can 
be decomposed into three arcs which cross the the centre 
of the cell at f(p), as shown in Fig. (J7J). Each of these 
three arcs can be further decomposed into three points 
which can be used to fit the equation of a circle, from 
which the curvature of the arc can be readily computed. 
Let us denote the three directly computed curvatures as 
Xi, Xi an d X3- We can define the mean square curvature 
as, 



1 



(xl + xl + xl)- 



(21) 



Thus, for the two conformal lattices shown in fig 4, we 
can estimate the mean square curvature for a given vertex 
using Eq. (f2"Tj) . The results are shown by the black dots 
in Fig. © and Fig. ©, along side are the theoretical 
results Eq. (19]) and Eq. ([20]) . 

Except in the case of complex inversion there will be a 
difference between the analytical value of the curvature 
and the estimated curvature. This is due to the fact that 
the curves to which we have fitted the circles are not 
themselves circular arcs. We quantify this discrepancy 
in section VI of this paper. 



8 




-0.5 0.5 

X coordinate 



U coordinate 




FIG. 6: Top left: a circular disk cut out of a triangular lattice in the z plane, where the angle between each successive heavy 
line is 7r/3. Clearly the lattice lines are free of curvature. Upon applying a conformal transformation, top right: w = z 1 ^ 2 
and bottom: complex inversion, i.e. w = 1/z, the result is a conformal lattice in the image plane. Note that in the case of 
w = z 1 ' 2 the conformal lattice has been produced by allowing the range of the z plane to extend to <f> < 4ir. In both cases 
the mean curvature vanishes everywhere but the mean square curvature does not. Also the mean square curvature diverges 
at the origin, for both transformations the origin is a critical point. At the critical point the angle preserving property of the 
transformation breaks down, this can be easily seen in the case of the w = z 1 ^ 2 transformation - instead of being preserved, 
the angle iv/3 between rays (heavy lines) emanating from the origin is halved. 
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FIG. 7: This diagram shows the effect of the transformation 
f(z) on a hexagonal lattice cell, which in the z plane has its 
centre at p. In the w-plane the lattice cell becomes deformed 
and its centre is now to be found at f(p). The deformed 
hexagonal cell can be decomposed into three arcs which cross 
the centre of the cell. Each arc can in turn be decomposed into 
points, to which we fit the equation of a circle. From this it 
is possible to calculate the curvature of each arc, represented 
here by xi,X2 and xs- 




FIG. 9: Mean square curvature for the mapping f(z) = 
1/z: this graph shows a comparison between the theoretical 
mean square curvature (from Eq. (20}), shown by the contin- 
uous line, and the mean square curvature computed directly 
from the conformal lattice using Eq. (|21[1 . shown by the cir- 
cles. 
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FIG. 8: Mean square curvature for the mapping f(z) = 
z 1 ^ 2 : this graph shows a comparison between the theoretical 
mean square curvature (from Eq. I19[0 . shown by the contin- 
uous line, and the mean square curvature estimated directly 
from the conformal lattice using Eq. (|2ip . shown by the cir- 
cles. 



V. SOME GENERALISATIONS 

In the previous section we considered a vertex, in the 
z-plane, to which there are connected a finite number of 
straight edges. We found that when the vertex is mapped 
to the image plane the mean curvature of the edges van- 



ishes and that the mean square curvature is independent 
with respect to the orientation of the vertex. 

Let us now consider the case where the edges connected 
to a vertex, in the z-plane, have some initial curvature. 
Let us also assume that the number of edges connected 
to the vertex is infinite. As the number of edges tends to 
infinity, the tangent vectors together describe a circle of 
unit radius centred on the vertex. Upon averaging over 
the length of the circle the mean curvature of the edges 
in the z-plane can be written as, 



C = 



1 

2^ 



(22) 



where k(4>) is now a continuous function over the range 
< (j) < 27r and describes the curvature of the edges in 
the z-plane. Similarly, the mean curvature of the edges 
upon being transformed to the w-plane is given by, 



C 



1 

27 



where we define, 



-Im 



hip) 



\h{p)\ L/i(p) 

Eq. (|2"4"]) can be written as 
1 



k(4>) 



\h(p)\ 
I/ 2 (p)I 
\h(p)\'< 



Im 



Im 



.,'<-> 



l/i(p)l' 



\hip)V 



(23) 



(24) 



l/i (p)1 
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where T and have the same definitions as in section 
IV. Writing out explicitly and simplifying gives 



K&) 



l/a(p)l 



sin(0 + 0) 



k{4>) 



I/iCp)I 2 v " ' r ' ' \fi(p)\' 

Substituting Eq. ([25]) into Eq. ((53J) yields, 

i |/ 2 (p)| ^ . _ . . f 



(25) 



27T |/l(p)| 2 7 

c 



sin(9 + - 



\fi(p)\ 



(26) 



Thus the conformal transformation scales the original 
mean curvature by the magnification factor l/|/i(p)| of 
the transformation (see Eq. (JTJ). 

Since the mean curvature does not necessarily vanish, 
we define the mean square curvature in the z-plane as, 



1 



2ji 



(27) 



similarly we define the mean square curvature in the w- 
plane as, 



Q = 



l 

2^ 



(«(0) 



C 



Expanding Eq. (|28|) yields, 



1 

2^ 



- 2Ck(0) + C J d<£, 



(28) 



(29) 



substituting Eq. (f25[) into Eq. (|29|) and performing the 
resulting integrals gives, 



Q = 2 



1 / I/ 2 (P)| 



Q 



I/i(p)I 2 
i I/ 2 (p)| J 



^i/i(p)i 3 y 



i/iwi 2 ' 

k(4>) sin(6 + 



(30) 



The first term is the mean square curvature of the edges 
in the image plane if the edges in the z-plane have no cur- 
vature. The second term simply states that the original 
mean square curvature of the edges is scaled by a factor 
of l|/i(p)| 2 . The third term is a coupling between the 
original curvature of the edges and the transformation 
/(*)• 

To get a better understanding of the third term let us 
label it as Q3. It can be written as, 



1 \h(p)\ 
^I/i(p)I 3 io 



2 71 



k (4>) (sin 9 cos (f) + cos sin 



(31) 

Since k(4>) is an arbitrary function which exists over the 
interval < <p < 27r, it is natural to express it in terms 
of a Fourier series, so that, 



1 



> a ° + a i cos (#) + b q sin(g0). (32) 
9 =i 



(33) 



Substituting Eq. (@2j) into Eq. (J3TJ) we find, 
03 = JT7% («i sin + fo i cos Q ) . 

\h(p)r 

thus the only contribution to the mean square curvature 
is from the mode 9 = 1. The contributions from all other 
modes cancel out, so that an increase in mean square 
curvature for a given edge is cancelled out by a decrease in 
mean square curvature for some other edge. Furthermore 
Eq. (|33[) depends on the orientation of the vertex in the 
z-plane, through the angle which is the argument of the 
complex function f2(p) / fi(p) and where p is the position 
of the vertex in the z-plane. The turning points (i.e. 
the location of the maximum and minimum mean square 
curvature) depend on the values of the co-efficients a± 
and 61 and are located at = tan~ 1 {ai/bi) and = 
tan~ 1 (ai/bi) + w. 

Thus when the edges in the z-plane are not free of cur- 
vature the effect of the conformal mapping is to scale 
the mean curvature, while the mean square curvature is 
found to be no longer independent with respect to orien- 
tation of the vertex. 



VI. HIGHER ORDER TERMS IN THE 
COMPLEX CURVATURE 

In this section we shall distinguish between the instan- 
taneous complex curvature of the image curve, which is 
given by Eq. ([3]), and the average complex curvature of 
the image curve - which we define below. 

For a given 2D plane curve its curvature is defined 
as the rate of change of the angle of the tangent vec- 
tor with respect to the distance along the curve. We 
can think of this as limiting process: take two tangent 
vectors a short distance apart on the curve, measure the 
angle they make with respect to a fixed axis, and take the 
limit in which the separation between the tangent vectors 
vanishes. This then gives the instantaneous value of the 
curvature. It is in this manner that the expression for 
the instantaneous complex curvature, given by Eq. (J5}, 
is derived. 

Consider a 2D plane curve in the image plane which has 
been generated by applying a complex mapping to some 
curve in the z plane. We can choose a pair of points on the 
image curve a finite distance apart and for both points 
construct a tangent vector. Thus we can ask: what is the 
rule which gives the total difference in angle between the 
two tangent vectors? 

In this section we derive an expression for the average 
complex curvature which we define as: the total change in 
the angle of the tangent vector over a segment of a curve 
averaged over the length of the segment. The result is 
a series expansion in which Eq. ([3]) is the lowest order 
term. This series expansion will lead us to an expres- 
sion which can be used to measure the degree to which 
a given image arc differs from being perfectly circular. 
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FIG. 10: Top: This diagram shows a straight line K in the 
z-plane. The line starts at the point labelled p and ends at 
q, it has a total length L. At both points the unit tangent 
vector has the same direction and makes an angle of <j> with 
the x-axis. Bottom: This diagram shows the resulting curve 
K in the image plane due to the conformal mapping w — f(z). 
The curve begins at the point labelled f(p) and ends at f(q) 
and now has a length I. In general the angle that the tangent 
vectors make at f(p) and f(q), which are labelled cj>\ and <f>2 
respectively, will no longer be equal. 



This is an important consideration for the applicability 
of conformal transformations in the context of foams. 

Here we shall only consider the effect of the mapping 
w = f(z) on a straight line in the z-plane. Furthermore, 
we shall assume that the effect of the mapping is to pro- 
duce a open, simple and well behaved curve in the image 
plane. 

Consider a straight line K in the z plane, which joins 



point p to q and let us assume it has a length |£|. This is 
shown in the top part of Fig. (|10[) . At p we have drawn 
the unit tangent vector £ = e 1 ^, so that the vector joining 
p to q is given by 

€ = (34) 

The line K is described by the equation 

K = z(t) = p + t£ = re w + te** ', (35) 

where < t < Note that at t — we have z = p and 
at t = |£| we have 

z(\H\) = q = P + ti = re* e + \li\e^. 

Thus K has been parametrised by the free parameter t 
which measures the distance along K from the point p. 
The line K has a real component 

x(t) = r cos 9 + t cos cj>, 

and an imaginary component 

y(t) = r sin 9 + t sin <p. 

We notice that since K is a straight line then the tangent 
vector is the same at both points p and q. 

Now consider what happens if we apply an analytical 
mapping w — f(z) to K, the result is a new curve K in 
the image plane which starts at the point labelled f(p) 
and ends at f(q) - as shown in the bottom part of Fig. 
(fTOjl . This curve now has some new length which we call 
I (not to be confused with the chord distance between 
f(p) and /(<?)). 

We define the average complex curvature of the image 
curve K as 

~ fa - <h 

K= — = — 1— (36) 

where the tangent vectors at f(p) and f(q) are labelled 
£i and £2 and make an angle of </>i and 02, respectively, 
with the the u axis. 

If for Eq. (|36[) we keep the point p fixed and take the 
limit in which |£| goes to zero we recover the instanta- 
neous curvature of the curve K at the point f(p), that 
is 

k(v) = Urn K, 

(note that «(p) refers to the instantaneous curvature of 
the image curve at the point f{p), in keeping with the 
definition of Eq. ([3]) we refer to the point f(p) by its 
corresponding z-plane coordinate; of course we are free 
to express k(j?) in terms of the image plane coordinates 
f(p) by making a change of variable - as we did for the 
example given in section II. B). 



12 



Let us now define the difference between the instanta- 
neous curvature of the image arc at its starting point and 
the average curvature of the arc, which we denote as 

AK = K-«(p). 

We note that if the image arc is a perfect circle then the 
instantaneous curvature has a constant value at every 
point on the arc and therefore |AK| = 0. This fact can 
be used to quantify the degree to which the image curve 
deviates from perfect circularity. 

In the context of 2D and quasi 2D foams it is impor- 
tant to know the degree to which a given arc, generated 
by a conformal mapping, deviates from perfect circular- 
ity. Laplace's law requires that the curvature of a soap 
film is constant, therefore for a 2D soap film a large value 
of |AK| implies a large error in how well the soap film 
approximates the conformal lattice. The same holds true 
for quasi 2D foams (in the constant volume regime), ex- 
cept that we have to remember that the total curvature of 
the soap film, i.e. the sum of longitudinal (n = /«(p))and 
transverse (at) curvatures, has an instantaneous value of 



H 



In either case Laplace's law is satisfied only if the soap 
film has a constant longitudinal curvature, i.e. it de- 
scribes a circular arc when viewed from a direction per- 
pendicular to the bottom plate of the Hele-Shaw cell. 

In the following we shall decompose the task of finding 
K into two parts: firstly we compute I and secondly we 
compute A(f>. 



A. Arc Length in the Image Plane 

As stated above: the real and imaginary components 
of the line K are a function of the free parameter t. Upon 
mapping K to the image plane, (using the analytical func- 
tion w = f(z)) the result is the image curve K - the real 
and imaginary components of which are also functions 
of the free parameter t and denoted by u{t) and v{t), 
respectively. We can write the real and imaginary com- 
ponents of K in the form of a series expansion about the 
point p (i.e. t = 0), giving 

u(t) = u + t^+t 2 ^L + ...+h.o.t. (37) 



and also 



where, 



v(t) = v + t^+t 2 ^ + ... + h.o.t. (38) 



u(t)\ 



t=0 



du{t) 
dt 



u(0) 



(39) 



(40) 



d 2 u{t) 



dt 2 



(41) 



t=o 



and so on. We also define v t and v tt (and higher order 
terms) in a similar manner. 

The total length of the arc K is given by, 




(42) 



We can substitute Eq.(37|) and Eq.(38J) into Eq. (42|) . then 
expand the integrand in powers of t, and integrate each 
term with respect to t to give, 



^o|£l+^| 2 + ^| 3 



h.o.t. 



(43) 



where we assume that |£| is small enough to guarantee 
convergence of the series. The coefficients in Eq. (|4"3"|) are 
defined as 



k = VV + 



vr 



1 U t Utt + V t Vtt 



1 



6 

{u t u tt + v t v a ) 



2 \Ju t 2 + v t 2 



Uu + UtUttt + VtVttt + vu 



(44) 



(45) 



Vu7 



2{ut 2 + v t 2 f 



(46) 



We notice that the coefficients given by Eq. (44|) to Eq. I46j) 
(and higher order terms) are defined with respect to the 
parameter t. In order to express the coefficients in terms 
of the complex number z we use Eq. ((35]) from which we 
have the relationship dz/dt = e^, this can be rearranged 
to give, 



dt = e'^dz. 



(47) 



Thus Eq.( [55|) to Eq.flHJ) (and higher order derivatives) 
can be written as 



d n u 



t=o 



Re 



Re 



Re 



d n w 



dt 



t=o 
d n w 



(e-^dz) 11 
d n w 



dz 11 
\fn(p)\Re 



in(j> 



,i(6 n +n<t 



|/„(p)| cos(0 n + n<j>) 



(48) 



where in going from the first line to the second line we 
have made a change of variable from t to z and used the 
fact that at t = we have z = p; also we used 



dz r 



= \fn(z)\e lAr ^^ = \.f n (p)\ 



z—p 
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Similarly for the derivatives Vt and vtt (and higher order 
terms), we have, 



h = ^Re 



hip) 



(51) 



IF 



\f„(p)\ sin(6>„ + ruj)). 



(49) 



Finally, the coefficients given by Eq.(|44p to Eq.([46p can 
be expressed in terms of Eq. (|48| and Eq. (|4"9")l to give 



*o = |/i(p)| 



(50) 



\fM\ \hi P )f 
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\h(p)? 



Re 



hip) 
MP) 



r 



ZDs 

hip) 



(52) 



Upon substituting Eq.®, Eq.® and Eq.® into 
Eq. (j43[) we have a series expansion in powers of the |£|. 
We note that the lowest order term in Eq. (4"3")l is identical 
to Eq.Q. 



B. Change in the angle of the tangent vector in the 
image plane 

We now compute the difference in angle, A<j> = 4>2~ </>i, 
of the tangent vectors in the image plane. This derivation 
is adapted from the one used by Needham to compute the 
instantaneous curvature (2j. The difference here is that 
we endeavour to retain all details which may yield higher 
order terms in the expression for A</>. 

Consider again the the straight line K in the z-plane, 
as shown Fig. (fit)]) . Note that the tangent vector at p 
upon being mapped by the function w = f{z) is, after be- 
ing translated to f(p), magnified by a length |/i(p)| and 
rotated by an angle Arg[fi(p)]. Similarly, upon applying 
the mapping w — f(z), the tangent vector at q is, after 
being translated to f(q), magnified by a length |/i(<7)| 
and rotated by an angle Arg[fx(q)]- However, the rota- 
tion of the tangent vector at f(q) will differ very slightly 
by A</> from the rotation suffered by the tangent vector 
at f(p). 

To compute the difference in angle between the tan- 
gent vectors consider the two points p and q as shown on 
the left hand side of Fig. (JTTJ) , upon applying the map- 
ping f\(z) the points are now located at hip) an d fiil), 
respectively. Where the vector connecting /i (p) to fi (q) 
is labelled x an d is given by 



x = hiq)-hip) 

= hi P )t+lihiP)e + lihiP)e- 



h.o.t 



hipMl 



^hip)^\e+^hip)^\e 



(53) 



+ ... + h.o.t 
where we have made use of Eq. (|3~4"|) . 

Note, with reference to the right hand side of Fig. (fTT)l . 
that to transform the vector fi(p) into the vector f\{q) 



it is necessary to expand the length of hip) by some 
factor and to rotate the vector fi(p) by an angle a. It is 
this angle a which is the extra rotation suffered by the 
tangent vector at f(q) compared to the rotation of the 
tangent vector at f(p), i.e. 

a = Acf. 

Our task is to compute the angle a shown on the right 
hand side of Fig. (fTTj) . Note, a is part of a triangle whose 
sides are given by the known vectors %, hip) and fi(q)- 
To compute a we rotate the triangle shown on the right 
hand side of Fig. (fTl~j) to the real axis by dividing by fi (p) 
(this is the same trick used by Needham Q), the result 
is shown in Fig. (|12[) , this new triangle is geometrically 
similar to the original triangle shown on the right hand 
side of Fig. (fTl~j) . The side along the real axis has a length 
1 and in addition the point v is given by the vector 

»=\x/hip)]- (54) 
By a simple application of trigonometry we have 

Im\v\ 



A(f> 



tan 



1 + Re M 



(55) 



Upon substituting Eq. (154)1 into Eq. (155)1 . we can expand 
Eq.([55]) in powers of |£| (we assume that |£| is small 
enough to allow convergence) to give, 

A4> = a = a \i\ + aiiei 2 + <t 2 |£| 3 ... + h.o.t. (56) 
where we have made use of Eq. ((53")) , we find 

'hip), 



Co = Im 



flip) 



-Im 



hip) (hip) 



hip) \hip) 



(57) 



(58) 
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FIG. 11: Left: This diagram shows two points p and q in the z plane. Right: Upon applying the mapping fi(z) — df(z)/dz 
the points are mapped to fi(p) and fi{q), which are connected by the vector \- Note that to transform fi(p) into fi(q) it is 
necessary to magnify /i (p) by some factor and to rotate fi (p) by the angle denoted by a 



<J<2 = Im 



' hip) 
JAP) 



(h(p)) 2 



hip) 
hip) 



W) 3 



(59) 




Im[v] 



FIG. 12: This triangle is the result of dividing all the vectors 
shown on the right hand side of Fig. JTTJ by the complex 
number flip)- Dividing by fi(p) has two effects: firstly the 
triangle is rotated so that one of its sides is now parallel to the 
real axis, secondly all the sides of the triangle are uniformly 
scaled. This second effect means that the new triangle is 
similar to the triangle shown in the right hand side of of Fig. 

Gil 



C. The Average Curvature 

Substituting Eq.([43]) and Eq.([56|) into Eq. flU and 
expanding in powers of |£| yields, 



where we find, 



\Q — K — 



\hip)\ 



111! 



hip) 
HP) 



.h.o.t 



(60) 



2\h(p)\ 



-Im 



Im 



hi?) 
MP) 



hip) 
Mp) 



i 



^2 = 



3|/ 2 (P)| ; 
*I/i(p)I V \hip)\ 2 



Im 



hip) 
hip) 



27Im 



hip) 
hip) 



lip)? 



\hip)\ 2 
341m 



Im 



Mp) 



8Im 



hip) 
hip) hip) gTTN 3 

-rmr m 

hip) 



In effect Eq. ^rjO]) gives the total change in the angle of 
the tangent vector in the image plane, as it traverses the 
image curve, averaged over the length of the curve. We 
observe that as |£| — * then K — > Ko = k, giving the 
instantaneous curvature of the image curve at the point 
/(p). In order to use Eq. lfiTJ)) we must define the following: 
(i) a straight line in the z-plane of length |£|, for which 



£ = e 1 * is the tangent vector at the point p 



re 



w. 



(ii) 
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FIG. 13: A straight line in the z-plane, with r = 1, = 0, (f> = 
tt/2 and |£| = 0.1 




FIG. 14: The resulting curve in the image plane generated 
by applying the conformal mapping f(z) = z 3 to the straight 
line defined by Eg. (161) 



by Eq.([35|). and let us set r 
|£| = 0.1, so that we have 



K = z(t) = l+te^, 



0,0 = tt/2 and 



(61) 



and < t < 0.1, the result is the straight line shown in 
Fig. (|13p . Upon applying the conformal transformation 
f(z) = z 3 the result is the curve shown in Fig. (fT4|) . 



The derivatives for the function f(z) 



- -,3 



are given 

by AO) = 3z 2 , f 2 {z) = 6z, f 3 (z) = 6, with all higher 
derivatives being equal to zero. Upon evaluating these 
derivatives at the point p — re 10 — le lQ — 1 we have 
fi(p) = 3, f2{p) = 6 and fz(p) = 6; together with the 
angle of the tangent vector to the straight line at the 
point p (i.e. <fi — tt/2) and the length of the straight line 
in z-planc (|£| = 0.1), we find 



-Im 



G 



itt/2 



• sin7r/2 



and similarly Ki = 0, K 2 = -1/2.25, K 3 = and K 4 = 
0.281481; thus up to fourth order in |£| we have IK « 
0.66225037. 

This value can be compared to a numerically computed 
value. To do so, for the curve shown Fig. (fT4")) . we need 
to compute the length of the curve and the total change 
in the angle of the tangent vector. This can be done by 
taking a large number of equally spaced points on the 
curve. By drawing a straight line between any two adja- 
cent points we get a series of segments which approximate 
the curve, the approximation improves as more points are 
taken. It is then a simple matter for a computer to cal- 
culate the length of each segment and to find the total 
length of the curve by summing up the lengths of all the 
segments. To compute the total change in the angle of 
the tangent vector: we first compute the angle that the 
first segment makes with the u axis, then we compute the 
angle that the final segment makes with the u axis, by 
taking the difference between the two we find the change 
in the angle of the tangent vector between the start and 
the end of the curve. Upon taking 1x10 s such equally 
spaced points on the curve shown in Fig. (|14[) we have 
an numerically computed value for the average curva- 
ture, giving K Nu m = 0.6622435. It is clear that Eq. ffeTJ]) 
is converging to Knuhi as more terms are taken. 



an analytical function f(z) which maps the line to the 
image plane (for which we can compute the derivatives 
f n (z) at the point p). 

Example f(z) = z 3 

Let us apply Eq. l60p to two simple examples, we choose 
the complex mapping f(z) = z 3 and inversion f(z) = 

Consider again the a straight line in the z-plane given 



Example f(z) — 1/z 

Upon carrying out the same procedure outlined above 
for inversion, i.e. f(z) — 1/z, also applied to the straight 
line Eq.f l6"Tj) (shown in the top part of Fig. (H31))- 
we find^IKo j= k = 2 and that all other coefficients 
Ki = K2 = K3...etc = 0. This demonstrates that for a 
curve produced by complex inversion the instantaneous 
curvature, at the starting point f(p) of the curve, is equal 
to its average curvature. This can only be the case if the 
curve has a constant curvature and is therefore a circular 
arc. 
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D. Mean Square of the Average Curvature 

For the instantaneous curvature we found it useful to 
compute the mean square curvature. The mean square 
curvature is independent of the orientation of the vertex 
and therefore serves as a useful measure of how strongly 
the conformal lattice is curved at any given point in the 
image plane. 

Similarly we can compute a mean square value of the 
average curvature. We find that this leads to a series 
expansion in which the lowest order term is the mean 
square curvature given by Eq. (I17|). When the higher 
order terms in this expansion vanish it means that all 
the image edges, at that particular point in the image 
plane, have a constant curvature. If on the other hand 
the higher order terms do not vanish this means that 
Laplace's law, for a 2D or quasi 2D foam used to realise 
the conformal lattice, cannot be perfectly satisfied for ev- 
ery edge. The magnitude of the higher order correction 
serve as an estimate of the discrepancy between the foam 
and the conformal map, in most practical cases the dis- 
crepancy can be assumed to be negligible if the leading 
order correction in the expansion is small. 

Consider again a vertex p in the z-plane to which there 
are connected a number of straight edges of length |£| 
separated by a constant angular separation, such as that 



shown in Fig. ((4]). Upon being mapped to the image 
plane the point p is now located at f(p) and the straight 
edges are transformed into arcs. If we were to sum up 
the average curvature, for each of the arcs emanating 
form the vertex at f(p), we would find that it vanishes - 
this is equivalent to saying that the mean of the average 
curvature always vanishes. 

If, however, for each arc we square the average curva- 
ture and sum up the squares we have a quantity that does 
not vanish. Thus we can compute the mean square value 
of the average curvature K, which we shall denote by Q, 
this can be done, most easily, by integrating Eq. (60]) with 
respect to <j> over the range < <j> < 2it giving 



where we find 



— / K(0) <#. 
27T Jo 



= Mo 



hici 2 + o(iei 4 ) 



I/2G01 



(62) 



i/i(p)p 



and 



J 



\f2( P )\ 4 
\fl(p)\ 2 \\Mp)\ 4 



32 12 



r 



/i(p)/ 3 (p) 



i\Mp)1 
8\MpW 



(63) 



There are two circumstances under which Q = Qo = Q- 
The first, which is the trivial case, is when the lattice 
spacing |£| of the original lattice in the z-plane vanishes. 
The second case is complex inversion f(z) = l/z (or more 
generally a bilinear transformation), for which we find 

Qi = Q2 = Qa--' = 0. The vanishing of all the higher or- 
der corrections means that all the curves emanating from 
the vertex at /(p) are circular arcs. It follows that, for 
a 2D foam or a quasi 2D foam, Laplace's law is satisfied 
by each of the arcs. 

Let us now turn our attention to the mapping f(z) = 
z 1 / 2 in this case we do not expect Laplace's law to hold 
for each of the arcs. We have already computed the low- 
est order term in Eq. ([62|) and found that 

Qo = Q = ^ 

(see Eq.QH)). Upon substituting f(z) = z 1 / 2 into Eq. f55|) 
we find 

~ _ J_l_ _ J_ 1_ 

< ^ 1 ~ 128^3 ~~ 128^6' 
where we have transformed from the coordinates in the 
z-plane (r, A) to the w-plane coordinates (R, A) by using 



the relationship *Jr = R. Assuming that the original 
lattice has a constant lattice spacing |£| we can say that 
a 2D foam, or a quasi 2D foam, is capable of producing 
a good approximation of a conformal lattice when 

- 

is small. We see that this happens at regions a large 
distance form central (critical) point of the conformal 
lattice, i.e when R » (7/128) 1 / e |CI 1/3 - 



VII. RELATED WORK 

Recently, the subject of conformal lattices with circular 
symmetry, of the type f(z) — z a , such as that shown in 
the top right hand part of Fig. ([6]), has been considered 
from the context of 2D crystals with non-uniform density 
by Mughal and Moore 

The authors performed simulations to find the ground 
state of a cluster of (classical) electrons confined to move 
on a disc. As expected, the charge density of the sys- 
tem is non-uniform, being small at the centre of the sys- 
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tern and very large towards the edge. The positions of 
the charges in the resulting ground state are triangulated 
using the Delaunay triangulation package Qhull, this al- 
lows the authors to uniquely define the number of nearest 
neighbours for each charge. The resulting triangulation 
is then visualised using the package Geomview, see for 
example Fig. (TT5]) - this image is from Q, 

By this method it is found that the ground state is 
crystalline and that most charges have six nearest neigh- 
bours. However, the cluster also contains a large number 
of defects in its interior. These defects are of two types: 
disclinations and dislocations. 

A disclination is any point which has an anomalous 
coordination number, that is any charge for which the 
number of nearest neighbours is not equal to six (or four 
on the cluster edge). A disclination is defined to have a 
+1 (-1) topological charge if its coordination number is 
equal to five (seven) - a similar rule applies to charges on 
the cluster edge. Topological constraints, due to Euler's 
theorem, ensure that the total topological charge in the 
system always sums to +6, i.e. there are always 6 more 
positive disclinations than negative ones. 

A dislocation, on the other hand, is defined as a tightly 
bound pair of disclinations of the opposite charge. Since 
dislocations are topologically neutral their numbers are 
not restricted by topology. Instead they are present in 
the lattice for energetic reasons. 

Disclinations induce an enormous elastic stress in the 
lattice. Consequently, in the ground state of 2D systems 
which form a triangular lattice with uniform density [l3| , 
the 6 positive disclinations demanded by topology are lo- 
cated on the edge of the lattice. In the case of systems 
with non-uniform density, Mughal and Moore [6( showed 
that disclinations are found in lattice interior. Specifi- 
cally, a system whose density increases (decreases) upon 
going from the centre towards the edge must contain 
more negative (positive) disclinations in its interior than 
positive (negative) ones. The total disclination charge 
must still sum to +6. In the case of the system shown in 
Fig. (fT5|) , it is found that for every negative disclination 
in the system interior their is a compensating positive 
disclination on the edge of the system. In addition, a 
further six positive disclinations are located on the edge 
to satisfy Euler's theorem. 

Mughal and Moore demonstrated that these interior 
disclinations are the mechanism by which the system 
maintains a crystalline structure while having a non- 
uniform density. However, the price that the system pays 
is the loss of translational and rotational symmetries, or 
in other words, the lattice acquires curvature. Using the 
framework of Burgers circuits and Burgers vectors, which 
is the natural way of characterising lattice defects, the 
authors are able relate the rate of change of the lattice 
density to the density of defects. In addition, the authors 
derive a relationship between the density of defects and 
the magnitude of the lattice curvature. 

During the course of this work the authors were led to 
consider conformal lattices. They suggested that the cen- 



tral point of the transformation of the type f(z) — z a , i.e. 
the critical point, should be considered as a disclination. 
The parallels between a disclination and the critical point 
in a conformal transformation of the type f{z) = z a are 
clear: both have an a anomalous coordination number, 
the lattice density is increasing/decreasing upon going 
away from the critical point/disclination and finally in 
both cases the surrounding lattice has a curvature. 

By treating the critical point, in a conformal lattice of 
the type f(z) = z a , as an isolated disclination at the cen- 
tre of the cluster, the authors were able to successfully 
compute the density of the conformal lattice (indepen- 
dently of Eq. ([T])) and the magnitude of the lattice cur- 
vature. It is satisfying to see that the expression given for 
the maximum lattice curvature in this earlier work agrees 
with Eq. (6|) . Conversely, the same theoretical arguments 
suggest that the ground state of a cluster containing in- 
ternal disclinations, such as that shown in Fig. (TT5]), is 
conformal everywhere except at the disclinations. 

A recent paper by Rivier et. al. Q also reaches 
the same conclusion: a disclination in a 2D crystalline 
structure can be regarded as a singularity, of the type 
f(z) = z a , in the complex plane. This connection be- 
tween conformal mapping and topological defects opens 
up an interesting avenue of research which we shall now 
discuss. 

Considerable effort has been devoted to finding the 
minimum energy configurations of interacting particles 
confined to 2D surfaces with a frozen topography. Sim- 
ulations and experiments have shown that the particles 
form an ordered triangular structure^in which most par- 
ticles have six nearest neighbours [l4j. In addition the 
ground state must also contain disclinations (and dislo- 
cations) . Once again the total disclination charge of the 
system depends on topology, for example: on a sphere 
the total disclination charge must sum to +12 while on 
a torus it is equal to (in these studies the choice of the 
interaction potential between particles leads to a uniform 
density of particles on the surface), see [l5[ and [l6l |. 
Such work has shown that even when disclinations are 
not required by topology, as in the case of the torus [l6[ , 
they can nonetheless be present in the ground state in 
response to the Gaussian curvature of the the underlying 
topography. 

The consideration of such problems, which are essen- 
tially concerned with finding the best way to uniformly 
arrange particles on a arbitrary surface, finds uses in the 
generation of computational meshes. The sites at which 
the particles are located, in the minimum energy config- 
uration, can be used to define a computational mesh for 
numerically solving partial differential equations. This 
technique has been used in the computational modelling 
of a host of physical phenomena, see for example [17| and 
references therein. Critically, the quality of the mesh for 
computational purposes depends on factors such as the 
mesh shape, size and smoothness. However, even for a 
modest number of particles a great deal of computational 
effort is required to find the optimal arrangement. 
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FIG. 15: From [g], with kind permission 5000 (classical) charges in a hard wall confining potential. The ground state 
of the system is crystalline with a non-uniform density. The density being smallest at the centre of the system and greatest 
towards the edge. The cluster interior contains disclinations, these are points whose coordination number is not equal to 6. 
A disclination has a positive (negative) topological charge if its coordination number is less (greater) than 6. In this image 
positive (negative) disclinations are coloured white (black). Note, as a consequence of the non-uniform density, the system 
contains more negative disclinations in its interior than positive ones. This excess of internal disclinations gives rise to the 
lattice curvature, which can be clearly seen towards the edge of the system. 



In response to this, more recent developments have 
shifted the focus from finding the optimal arrangement 
of the particles to finding the optimal arrangement of the 
defects themselves. This approach, pioneered by Bowick 
et. al. [18], integrates out the underlying triangular lat- 
tice to directly make the defects the degrees of freedom. 
The result is an elastic Hamiltonian, which can be numer- 
ically optimised to give the optimal configuration of the 
defects. This approach has the advantage of reducing the 
number of degrees of freedom and has been successfully 
used to find the arrangement of defects on the surface of 
a sphere [Hj]. The price that we pay for this decrease in 
computational burden is that the details describing the 
underlying lattice are lost. 

Knowing the position of the defects in itself is not suf- 
ficient information to use this algorithm for generating 
computational meshes. This is where conformal lattices 



might be useful. The idea is that once the optimal num- 
ber and arrangement of disclinations in the system has 
been determined using the algorithm outlined in [181 ] , the 
underlying lattice can be reconstructed by using a series 
of conformal transformations of an initial regular lattice. 
This approach may eventually lead to a more intelligent 
design of computational meshes. This is a topic currently 
being investigated by the present authors. 

VIII. CONCLUSION 

As detailed in Q , the conformal transformation of the 
honeycomb structure can produce ideal, dry, foam struc- 
tures that approximate some that can be realised in the 
laboratory. But the true equilibrium structure is slightly 
different, in that the local curvature is required to be 



19 



constant along each edge. This discrepancy sets limits 
on the applicability of conformal transformations in the 
context of foams. The results presented in section VI 
will be useful to experimentalists working in this field to 
estimate the magnitude of this discrepancy. 

Higher order terms for all the expansions 
given in this paper are available from AM 
(adil.m.mughal@gmail.com) . 
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